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Abstract 

Exact mean field equations are derived analytically to give the fundamental 
diagrams, i.e., the average speed - car density relations, for the Fukui-Ishibashi 
one-dimensional traffic flow cellular automaton model of high speed vehicles 
(vmax = M > 1) with stochastic delay. Starting with the basic equation de- 
scribing the time evolution of the number of empty sites in front of each car, 
the concepts of inter-car spacings longer and shorter than M are introduced. 
The probabilities of having long and short spacings on the road are calculated. 
For high car densities (p > 1/M), it is shown that inter-car spacings longer 
than M will be shortened as the traffic flow evolves in time, and any initial 
configurations approach a steady state in which all the inter-car spacings are 
of the short type. Similarly for low car densities (p < 1/M), it can be shown 
that traffic flow approaches an asymptotic steady state in which all the inter- 
car spacings are longer than M — 2. The average traffic speed is then obtained 
analytically as a function of car density in the asymptotic steady state. The 
fundamental diagram so obtained is in excellent agreement with simulation 
data. 

PACS numbers: 64.60.Ak, 05.40.+j, 05.70.Jk, 89.40.+k 
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I. INTRODUCTION 



Recently, there has been much interest in studying traffic flow problems within the con- 
text of cellular automaton (CA) models [1-4]. Compared with the fluid dynamical approaches 
to traffic flow problems, the CA models are conceptually simpler, and can be readily imple- 
mented on computers. These models capture the complexity of the nonlinear character of 
the problem and provide clear physical pictures [2,5-7]. For example, CA models show the 
existence of a transition between a moving phase and a jamming phase in the traffic of a 
city as the car density is varied [4]. These models also have the advantages that they can be 
easily modified to deal with the effects of different kinds of realistic conditions, such as road 
blocks and hindrances, traffic accidents [8], highway junctions [9], overpasses [10], vehicle 
acceleration [11], quenched disorderness [12], stochastic delay due to driver's reactions [13], 
anisotropy of car distributions in different driving directions [14], faulty traffic lights [15], 
etc. Recently, CA models have been successfully applied to study traffic flow in a city by 
performing high speed simulations on the actual road map of the city of Dallas [16]. In 
view of the increasing importance of CA models in studying traffic flow problems, it is thus 
important to understand these models in more detail, especially from the point of view of 
statistical mechanics and nonlinear dynamics. 

The basic one dimensional (ID) CA model for highway traffic flow [3,13] is the CA rule 
184 [1]. This model describes single- lane traffic on a road of length L with periodic boundary 
condition. Each of the L sites can either be empty or occupied by one vehicle. Let N be 
total number of cars, then the average vehicle density on the road is p = N/L. The cars 
move from the left to the right according to the following rules. All the cars attempt a move 
of one step to the right simultaneously at each time step. If the site in front of a car is not 
occupied at that time step, the car moves one site ahead. Otherwise, the car cannot move. 
This simple model predicts a transition from laminar traffic flow to start-stop wave as the 
car density increases. 

There are many variations on the basic model. Nagel and Schreckenberg (NS) considered 
the effects of acceleration and stochastic delay of vehicles with high speed [3,13]. In the NS 
model, a car can move at most by M sites in a time step. The actual speed at a time step 
depends on the spacing in front. If the speed in the present time step is less M and the 
spacing ahead allows, then the speed increases by one unit in the next time step. If the 
spacing ahead is less than the speed in the present time step, then the speed is reduced to 
the value allowed by the spacing, and thus leads to a deceleration. In addition, there is a 
probability that the speed of a car is reduced by one unit in the next time step. Thus the 
NS model captures the features of gradual acceleration, deceleration and randomization in 
realistic traffic flows. 

Fukui and Ishibashi (FI) introduced another variation on the basic model [17] in which 
the cars can move by at most M sites in one time step if they are not blocked by cars in 
front. More precisely, if the number of empty sites C in front of a car is larger than M at 
time t, then it can move forward M (M — 1) sites in the next time step with probability 
1 — / (/). Here, the probability / represents the degree of stochastic delay. The / = 
model is referred to as the deterministic FI model, while the / = 1 case is the deterministic 
FI model with the maximum speed v max = M — 1. If C < M at time t, then the car can 
only move by C sites in the next time step. The FI model differs from the NS model in 
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that the increase in speed may not be gradual and stochastic delay only applies to the high 
speed cars. Obviously, the two models are identical for M — 1. 

Fukui and Ishibashi have performed numerical simulations on the model [17]. The fo- 
cus in traffic flow problems is the so called fundamental diagram, which is the relationship 
between the average speed in the steady state and the car density. While numerical simula- 
tions provide us with accurate fundamenetal diagrams, it is desirable and useful to have a 
better qualitative understanding of the numerical results within some analytic approaches 
such as mean field theories. Numerical studies and mean field theories have been extremely 
useful in providing detailed understanding in phase transitions and critical phenomena in 
equilibrium statistical mechanics, and we foresee that they will be equally useful in the study 
of dynamical systems such as the CA traffic flow models. 

Various mean field theories have been proposed for traffic flow models in 2D [10,14,19-21] 
and ID [13,17,22-30]. In ID, mean field approaches giving results in exact agreement with 
simulations have been given for the NS model [13,22] with M = 1 and for the deterministic 
FI model [17]. In 2D, mean field theories have been proposed [10,14,20,21] for the model 
introduced by Biham et al [4] (BML model). Most of the mean field approahces are macro- 
scopic theories in that the consideration is based on the idea that the average duration that 
a car stays on a site, while depending on the speed of the cars, determines the blockage on 
the car behind it and thus, in turn, determines the average speed. 

Recently, the steady state of CA traffic flow models have been studied within statisti- 
cal mechanical approaches [23-27]. While these studies are also mean field in nature, the 
approaches are based on microscopic consideration focusing on the time evolution of the 
occupancy on each site of the road. A nonlinear mapping between the macroscopic average 
speeds at two consecutive time steps can then be derived by carrying out suitable statisti- 
cal averages on the microscopic relations. The stable fixed point of the mapping gives the 
steady state average speed as a function of car denisty. For the deterministic Fl model, 
results in exact agreement with numerical data have been obtained [26,27]; while for the FI 
model with delay, the microscopic approach gives results in good agreement with simula- 
tions. Such microscopic theory has the advantage that it provides a systematic approach for 
the derivation of mean field results for the steady state. 

An alternative microscopic approach based on the time evolution of inter-car spacings 
have also been proposed recently for the FI model [28-30]. The idea is similar to the car- 
oriented mean field theory (COMF) [22]. In Ref.[30], we studied the deterministic FI model 
for arbitrary M and found that the inter-car spacings self-organize themselves into either 
long or short spacings in the steady state depending on the car density on the road. The 
fundamental diagram so obtained is in exact agreement with simuations. In this paper, 
we generalize this approach to study FI model with arbitrary M and arbitrary degree of 
stochastic delay /. We are able to derive a general expression for the average car speed in 
the steady state as a function of M and / which is in excellent agreement with numerical 
data. 

The plan of the paper is as follows. In Sec. II, we present the basic evolution equation for 
the inter-car spacings. The concepts of long and short spacings are introduced. The average 
speed is expressed in terms of the probabilities of finding long and short spacings. Section 
III gives a detailed discussion on how the inter-car spacings evolve in time according to the 
inter-car spacing of the car in front. Using the results in Sec. Ill, Sec. IV gives the steady 
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state result in the high car density regime. It is argued that all the spacings evolve into 
short spacings in the steady state. In Sec. V, we derive the average speed in the low density 
regime using the idea of detailed balance. Results are discussed in Sec. VI and compared 
with numerical data for M = 2 and M = 3 with arbitrary degree of stochastic delay. The 
extension of the present approach to other traffic flow problems is also discussed. 

II. THE DEPENDENCE OF AVERAGE TRAFFIC SPEED ON INTER-CAR 

SPACINGS 

Let C n (t) be the number of empty sites in front of the n-th car at time t. It is also 
the distance between the n-th car and the (n + l)-th car. The average distance between 
neighboring cars can be represented by C = ^j^- — ^ — 1. Let v n (t) be the number of sites 
that the n-th car moves during an update at the time t, i.e. the update between time t and 
t + 1. The number of empty sites in front of the n-th car at time t + 1 is 

Cjt + 1) = Cjf) + v n+l (t) - v n (t) (1) 

Within the generalized FI traffic flow model with the maximum car velocity v max = M 
and a stochastic delay probability /, the relationship between the velocity of the n-th car 
and the inter-car spacing ahead at time t is 

Vn (t) = F M (f,C n (t)), 

where 

\ c ' 

F M (f,C) = < M - 1, with probability / 
[ M, with probobility 1 

The average speed at time t is 

^)4e^,w4( e c n (t)+ e (M-f)\. (4) 

JV n=l JV [c n (t)<M-l C n (t)>M ) 

The sums in Eq.(4) correspond to sums over two different types of car-spacings. An inter- 
car spacing is labelled a long spacing if it consists of M or more sites, i.e. if C n (t) > 
M; while an inter-car spacing is labelled a short spacing if it consists of M — 1 or less 
sites, i.e. C n (t) < M — 1. Let N m (t) be the number of cars at time t with m empty 
sites ahead. The probability that a car is found to have a spacing of m sites ahead is 
give by P m (t) = N m (t)/N. The number of long spacings Ni ong (t) at time t is given by 
Niong {t) = J2m>M N m (t), and the number of short spacings N s h ort (t) at time t is given by 
N short if) = Hm=l N m (f). The probability Pi on gf) °f an inter-car spacing belonging to a long 
spacing is Pi ong (t) = Ni ong (t) /N; while and the probability P s hortf) of an inter-car spacing 
belonging to a short spacing is P s hort f) = N short (f) /N. Thus, the average car speed in 
Eq.(4) can be expressed in terms of P m (t) as 

M-l 

V(t)='£m- P m (t) + (M-f) P long {t) . (5) 

m=l 



(2) 



-/ 



if C < M - 1 
if C> M. 



(3) 
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III. TIME EVOLUTION OF INTER-CAR SPACINGS 

For the FI traffic flow model with stochastic delay, the inter-car spacings evolve in time 
in the following ways depending on whether the spacing is short or long and on the nature 
of the inter-car spacing of the car in front. Suppose the spacing of the n-th car is short, i.e. 
C n {t) < M - 1. If C n+1 (t) < M - 1, then 

C n (t + l) = C n+1 (t)<M-l. (6) 

If C n+ i(f) > M, then 

r a + -\ \ = \ M ~ l > with P robabilit y / ( 7 \ 

n{ + ' \M, with probability 1 - /. 1 ' 

Equation (6) follows from Eqs.(l)-(3) that for C n (t) < M — 1 and C n+1 (t) < M — 1, 
C n (t+1) = C n (t) + C n+ i(t) — C n (t) = C n+ i(t) < M — l. For Eq.(7), we introduce a stochastic 
Boolean variable 9 n (f) describing the stochastic delay of the n-th car with probability /: 

n i f\ — J I? with probability / , , 

n{}) = \ 0, with probability 1 - /. [ } 

It then follows from Eqs.(l)-(3) that for C n (t) < M - 1 and C n+1 (t) > M, 

C n (t + 1) = (M — l)9 n+1 (f) + M9 n+1 (l - /), (9) 

and hence Eq.(7). 

Suppose the spacing of the n-th car is long, i.e. C n (t) > M. If C n+1 {t) < M — 1, then 
Cn+iW < 6„(t + 1) < j Cn{t) _ x w . th probabUity 1 _ f (10) 
If C n+1 (t) > M, then 

!1, with probability f(l — f) 

0, with pobobility f 2 + (1 - /) 2 (11) 

— 1, with probability f(l — f). 

The proof of Eq.(10) goes as follows. It follows from Eqs.(l)-(3) that for C n (t) > M and 
C n+1 (t)<M-l, 

C n {t + 1) = C n (t) + C n+1 (t) - \{M - l)6 n (f) + M6 n {l - /)]. (12) 
Since C n (t)> M, C n (* + 1) > C n+1 (t). For C„ +1 (0 < M — 1, 



C n+1 (f) - [(M - l)0 n (/) + M0 n (l - /)] < 
Hence Eq.(12) implies 

C n (* + 1) < 



0, with probability / 

— 1, with probability 1 — /. 



C n (t) with probability / 

C n (t) — 1 with probability 1 — /, 
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and Eq.(lO) is proven. 

For Eq.(ll), from Eqs.(l)-(3) and the conditions C„(t)> M and C n+ i(t)> M, we have: 

C n (t + 1) = C n (t) + [(M - l)0 n+1 (/) + M9 n+1 (l - /)] - ([(M - l)0 n (/) + M0 n (l - /)]) 
= C n (t) + [M9 n+1 (l — f) — (M — l)9 n (f)\ + [(M - l)9 n+1 (f) — (M — l)9 n (f)\ 
+[M9 n (l - f) - M9 n (l - /)] + [(M - l)9 n+1 (f) - M9 n (l - /)] 

!1, with probability f(l — f) 

0, with probability f 2 + (1 - /) 2 ( 13 ) 

— 1, with probability /(l — /), 

which is Eq.(ll). 

Equations (6), (7), (10), and (11) give the time evolution of C n (t). Coupled with Eq.(4) 
or Eq.(5) gives the time evolution of the average speed V(t) and hence the asymptotic limit 
can be studied. 



IV. HIGH DENSITY CASE 

For high car densities (p > 1/M), the average inter-car spacings satisfies C = 1/ p — 1 < 
M — 1. It can be argued that in the asymptotic steady state of traffic flow, all inter-car 
spacings become short spacings, i.e., 

C n (t)<M- 1, Vn. (14) 

/,From Eq.(6), if every spacing is shorter than M, then 

C n (t + l) = C n+1 (t), Vn, (15) 

which implies that the spacing in front of the n-th car at t + 1 is simply the spacing of the 
n + 1-th car at time t. Thus as time increases, the traffic evolves as a continuous shift in the 
numbering of cars. Therefore, the situation corresponding to all the spacings are short is a 
steady state. Since 

L-N = J2C n (t) < (M-1)N, 

n 

the condition in Eq.(15) holds for p > 1/M. Under this condition, the average speed V is 
simply the number of empty sites divided by the total number of cars: 

v = = ~p~ 1 ' for p - l ^ M ' ^ 

The time evolution of C n (t) (Eqs.(6), (7), (10), (11)) ensures that for p > 1/M the 
steady state in which all inter-car spacings are short spacings is approached asymptotically 
regardless of the initial state of the traffic flow. The proof has been given [29,30] for the 
deterministic FI model with arbitrary v max = M. While a similar proof can be given, 
we simply note that the stochastic delay becomes ineffective at high car densities and the 
system behaves increasingly as a deterministic model in the high density regime. Moreover, 
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the proof in Refs. [29,30] works both for the deterministic model corresponding to / = and 
for the totally delayed model with / = 1 corresponding to a deterministic model with M — 1. 
The existence of stochastic delays will only affect the time for the traffic to approach the 
asymptotic limit from its initial configuration, but not the nature of the asymptotic steady 
state. 

V. LOW DENSITY CASE 

For low car densities (p < 1/M) and M > 2, the average inter-car spacing is C — 
l/p — 1 > M — 1. It can be shown that in the asymptotic steady state, every spacing will 
not be shorter than M — 1, i.e., 

C n (t)>M- 1, Vn (17) 

or equivalently, N — Ni — ... — N M ^ 2 = 0, where N m is the number of cars with m 
empty sites ahead. For the deterministic FI model with v max —M, it has been proven 30 that 
the steady state corresponding to C n (t) > M for all n in the low car density regime with 
p < v 1 +1 = j^q-j- is approached after a finite period of time. As both the / = and f = 1 
limits of the generalized FI model correspond to deterministic FI models with v max = M 
and v max = M — 1, respectively, the inequality in (17) holds. Thus, in the steady state, 

P = Pi = - = Pm-2 = 0. (18) 

It follows that Piong = 1 — Pm-i, and the average speed in the steady state (Eq.(5)) can be 
written as 

V = (M — 1)P M _! + (M — f)P long — (M — f) — (1 — f)P M -i. (19) 

Thus the problem of finding V amounts to obtaining Pm-i in the asymptotic limit. 

To obtain Pj in the steady state, we introduce Nj^ j± i to describe the number of inter-car 
spacings with their lengths changed from j at time t to j ± 1 at time t + 1. The probability 
of finding an inter-car spacing with length j at time t and length j ± 1 at time t + 1 is 

W^ j±1 (t) = Nj^ j±1 (t) /N = (Nj (t) /N) ■ [Nj^ j±1 {t)/Nj{t)\. (20) 

^From Eqs.(6) and (7), we have 

Nm-i^m (t) /N M -i (t) = (1 - f)(P M + Pm + i + ■■■), (21) 

and 

Wm-i^m (t) = (1 - /) P M -i • Piong- (22) 
Similarly, from Eqs.(10) and (11), we have 

Nm^m-i (t) /N M (t) = (1 - /) P M -i + (!-/)/ (Pm + P M+ i + ■■■), (23) 

and 
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(t) = (1 " /) P M (P M -1 + fPlong). (24) 

For Wj^ j±1 with j > M, Eq.(12) states that for C n {t) =j>M and C n+1 (t) = M - 1, 

C n (t + 1) = j + M - 1 - [(M - l)0 n (f) + M9 n (l - /)] 

= A(/) + (j-lR(l-/). (25) 

Similarly for C n (t) =j>M and C n+ i(*) > M, Eq.(13) gives 

C n (t + 1) = J + [(M - l)0 n+1 (/) + M0 n+1 (l - /)] - [(M - l)9 n (f) + M0 n (l - /)] 

= (j - i)o n (i - f)e n+1 (f) + (j + i)0 n (/)0„ +1 (i - /) 

+j[0 n (f)0 n+1 (f) + B (1 - f)e n+1 (l - /)]. (26) 

Hence 

Nj^M/Njit) = (1 - /) • [Pm-i + /(Pm + Pvm + ...)], (27) 

and 

Nj^tyNjit) = /(l - /) • (P M + Pm+i + -..)• (28) 
The probabilities Wj->j±i(t) are then given by 

W j ^ j - 1 (t) = (l-f)P j (P M - 1 + f-P long ), j>M (29) 

and 

W j ^ j+1 (t)=f(l-f)P r P lmg , j>M. (30) 
The asymptotic steady state should satisfy the detailed balance conditions given by 

W M -l^M(t) = W M ^ M -i(t)] 
W j ^ j - 1 (t) = W j - 1 ^(t), 3 >M. (31) 



Substituting Eqs.(22), (24), (29) and (30), we obtain 

Pm 1 — Pm-i 



(32) 



Pm-i / + (1-/)Pm-i' 
and 

■p- = p f{ lm M i i s J > M - (33) 

where a is a constant depending only on Pm-i- Equations (32) and (33) imply that 

P M = (a/f)P M - 1 , P M +i = aP M , -., P n = a n ' M P M - (34) 
Together with the equality 



1 oo 

C = --l= £ J'Pj, (35) 

P j=M-l 



we obtain 



c-f-(i- a y 

M-l 



(36) 



a 2 + M ■ a(l - a) + (M - 1) • /(l - a) 2 ' 
Substituting Eq.(33) for the constant a, we arrive at a quadratic equation 

(1 - f)P M _, + (C - M + 2/)P M _! -/ = (37) 
for Pm-i which gives a non-negative root 

M-C-2/ + v /(C'-M + 2/) 2 + 4/(l-/) 



M-l 



2(1-/) 



(38) 



Substituting Eq.(38) for Pm-i back into Eq.(19) for the average speed in the steady state, 
we finally obtain 



M-C-2/+ J(C -M + 2f) 2 + 4/(1 - /) 
V(t^ oo) = M - f ^ 



M-l + l/p- y/(l/p -1-M + 2f) 2 + 4/(1 - /) 
- g • (39) 

Equations (16) and (39) are the main results of the present work. They give the average 
speed of cars in the steady state over the whole range of car densities for arbitrary maximum 
velocity and degrees of stochastic delays. 

VI. DISCUSSION 

The fundamental diagram, i.e. the speed - car-density relation, of the FI model with 
stochastic delay is 

f M-l+l/p- A /(l/p-l-M+2/)2+4/(l-/) ! 
V(t -> OO) = I ] 2 — P — ~M (40) 

The general features of the speed in the steady state are that for given M and in low density 
regime (p < 1/M), different values of / correspond to separate curves; while in the high 
density regime (p > 1/M), the curves corresponding to different values of / coincide with 
each other and fall into one curve. The latter feature reflects that stochastic delays become 
ineffective for sufficiently large car densities and systems with different values of / behave 
in the same way. As the car density increases, the curves for different values of / meet at 
p = 1/M. 

In Ref.[17], numerical results have been reported for the M = 2 FI model without good 
analytic explanations. It is thus illustrative to compare the present result with numerical 
simulations. For M — 2, our general result (Eq.(40)) reads 
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In order to compare with the analytic result, we carried out numerical simulations on one- 
dimensional chain with 1000 cars and the length of the chain was adjusted so as to give 
the desired density of cars. Periodic boundary condition was imposed. The motion of the 
cars was followed and the average speed of cars recorded. The first 20000 time steps were 
not included into the averaging procedure so as to ensure that the system has approached 
the steady state. The averages were taken over the next 80000 time steps. Our numerical 
results are consistent with those reported in Ref.[17]. Figure 1 compares the analytic results 
with numerical results for different values of /. Excellent agreement is found. Equation (41) 
thus complements the numerical results in Ref.[17] and provides an analytic expression for 
the numerical data. 

To further establish the validity of our result, we carried out numerical simulations for the 
M = 3 FI model with stochastic delays and compared results with our general expression. 
Results are shown in Fig. 2. Again, it is obvious from the figure that the agreement is 
excellent. In passing, we note that the M = 1 results agree well with numerical data. For 
M — 1, the FI model is identical to the NS model with a maximum velocity of unity and the 
M = 1 results reported here agree with those reported in the literature for the NS model 
for this particular case. 

In summary, we derived exact results for the average speed in the asymptotic steady state 
as a function of the car density and the degree of stochastic delay for the Fukui-Ishibashi 
traffic flow model. The approach is based on the study of the time evolution of the inter-car 
spacings. The notions of long and short inter-car spacings are introduced. The probability 
of finding a long or short spacing on the road is then calculated. In the high density regime 
(p > 1/M), all inter-car spacings will become short spacings in the steady state; while in 
low density regime (p < 1/M), all inter-car spacings will be longer than or equal to length 
M-l. The probabilities that a spacing becomes longer and shorter by one unit of length in a 
time step are calculated. The asymptotic steady state is obtained by imposing the condition 
of detailed balance. The general expression for the average speed in the steady state is then 
obtained analytically as a function of the car density for arbitrary value of the maximum 
velocity and arbitrary degree of stochastic delay. Results are compared with numerical data 
for M = 2 and M = 3 over the whole range of < / < 1. Our analytic results are in 
excellent agreement with numerical data. 

The present approach provides an alternative way to study traffic flow problems analyt- 
ically. In principle, our approach can be extended to study other models in one and two 
dimensions [3,4,9,10,13,14,21]. While it is relatively simple to study the time evolution of 
car spacings in one-dimensional models, it is non-trivial to extend the present approach to 
two-dimensional models. In 2D model such as the BML model [4], the spacings along one 
direction will be coupled to the evolution of the spacing in another direction as a group of 
short spacings in one direction will slow down the traffic flow in the perpendicular direction 
and hence influence the car spacings in the other direction. This coupling leads to compli- 
cated coupled equations for the probabilities of having long and short spacings in the two 
directions. Although more complicated than the one-dimensional case, this coupling in turns 
leads to the interesting phenomena of having a jamming to moving phase transition at a 



10 



finite car density. Work along this line is in progress and results will be reported elsewhere. 
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Figure Captions 

Figure 1: The fundamental diagram of the Fukui-Ishibashi traffic flow model with the 
maximum car velocity M = 2 and for different values of the degree of stochastic delay /. 
The solid curves are the theoretical results in Eq.(41). The points with different symbols 
represent results obtained by numerical simulations. The curves from the top down along 
the average velocity axis correspond to different values of / between / = to / = 1 in step 
of 0.1. 

Figure 2: The fundamental diagram of the Fukui-Ishibashi traffic flow model with the 
maximum car velocity M = 3 and for different values of the degree of stochastic delay /. 
The solid curves are the theoretical results in Eq.(40). The points with different symbols 
represent results obtained by numerical simulations. The curves from the top down along 
the average velocity axis correspond to different values of / between / = and / = 1 in 
step of 0.1. 
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